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Abstract 

The central limit theorem states that the sum of independently distributed n- 
tuples of real variables (subject to appropriate normalization) tends to a multivariate 
gaussian distribution for large A^. Here we propose to invert this argument: given a set 
of n correlated gaussian variables, we try to infer information about the spectrum of 
the discrete microscopic probability distributions whose convolution generated such a 
macroscopic behavior. The techniques developed along the article are applied to prove 
that the classical description of certain macroscopic optical experiments is infinitely 
more complex than the quantum one. 

1 Introduction 

The central limit theorem is one of the most celebrated results in statistics and applied 
probability. Originally formulated by Laplace, it states that the sum of several independent 
variables tends to a gaussian distribution in the asymptotic limit. The central limit theorem 
is used to estimate error bars in all sorts of statistical studies. It also justifies why the 
probability density of magnitudes as diverse as the IQ of a large sample of individuals or 
the noise in a radio signal approaches the so-called 'Bell curve', by postulating that such 
quantities are the sum of many independent (microscopic) contributions. 

In this work, we carry this reasoning one step further and consider the problem of extract- 
ing information about such 'mysterious microscopic variables' whose convolution generated 
the experimental macroscopic data. More concretely, we are interested in how different 
constraints over discrete microscopic models manifest in the macroscopic limit. 

The reason why we restrict our study to discrete models is two-fold: on one hand, from 
a practical point of view, discrete systems adopting a finite (small) set of possible values 
are preferred in computational modeling, due to the modest memory resources required to 
store and manipulate them. If the internal mechanisms of a given -say- biological process at 
the cellular level are unknown but there is plenty of data available about the macroscopic 
behavior of a group of cells, it is thus advisable to search for discrete models for individual cell 
behavior. On the other hand, from a foundational point of view, there exist some physical 
quantities, like time and space, which, although traditionally regarded as continuous, could 
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actually be discrete. It is thus interesting to know if such a discretization leaves a signature 
at the macroscopic level. 

In this paper we will perform a thorough analysis of this mathematical problem. We 
will show that, if the spectrum of the microscopic variables is finite and fixed a priori, then 
there always exist macroscopic bivariate distributions which cannot be approximated with 
such microscopic models. In fact, deciding which gaussian distributions are approximable 
or not can be cast as a linear program. We will also arrive at the surprising conclusion 
that certain spectra of infinite cardinality do not allow to reproduce all bivariate covariance 
matrices. If the spectrum of the microscopic systems is free but finite, then several n-tuples 
of n + 1-valued variables are needed to recover all n-variate gaussian distributions, and the 
problem of deciding if a given gaussian distribution is generated by rf-valued microscopic 
systems can be solved via semidefinite programming. 

Based on the ubiquitous applicability of the central limit theorem, we expect our results to 
be applicable in a large variety of scientific situations where one wants to infer microscopic 
information from macroscopic data, such as molecular biology, genetics, neuroscience or 
social sciences. We will illustrate this potential in the context of quantum mechanics. 

In particular we will apply the former results to show that light is infinite dimensional, 
or formally, that homodyne measurements over bipartite gaussian states of light cannot be 
understood as collective spin measurements over ensembles of independent pairs of correlated 
microscopic particles. We will also shed new light in the Foundations of Quantum Mechan- 
ics: in [16] it was shown that any quantum experiment admits a classical explanation when 
brought to the macroscopic scale (this phenomenon was named macroscopic locality). How- 
ever, the particular classical model found in [16] seemed excessively convoluted: whereas 
the quantum variables could adopt a finite number of values, the spectrum of the classical 
ones was continuous, ranging from — oo to +00. Using the techniques developed along this 
paper, we will show that no finite classical model can account for the intensities observed 
when both sides of a bipartite state consisting of many copies of the maximally entangled 
state are subject to extensive equatorial spin measurements. Since each of the microscopic 
quantum variables involved can adopt one of two values (up or down) in each instance of the 
experiment, the quantum description can thus be considered infinitely less complex than its 
classical counterpart. 

The paper is structured as follows: first, we will describe in detail a simple scenario where 
microscopic effects over the macroscopic variables can already be seen. We will give a first 
illustrative application of that to Brownian motion. Then, in Section [31 we will study the 
general problem of characterizing gaussian distributions generated by microscopic variables 
with a fixed output structure. We will treat both the finite and infinite dimensional case. 
We will also address some problems which may arise in practical implementations of macro- 
scopic experiments. In Section H] we will analyze the case where the set of possible values of 
the microscopic variables is not known a priori and may actually vary between independent 
tuples. The complexity of the problem will be treated next, and the consequences for Foun- 
dations of Quantum Mechanics will be covered in Section [61 In Section [7] we will present our 
conclusions. 
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2 A simple example 



Suppose that two separated parties, say, Alice and Bob, perform an experiment like that of 
Fig. [H that allows them to measure the continuous variables X and Y, respectively. After 
many repetitions of the experiment, they find that X and Y follow a bivariate gaussian 
probability distribution, defined by the values 

(X), {Y)nxx ^ {X') - (X)^7ry = {Y') - {Y)\ixy ^ {XY) - {X){Y). (1) 

Alice and Bob then make the hypothesis that X and Y are not gaussian by accident, but 
because they actually correspond to the sum of many pairwise-correlated independent vari- 
ables. A plausible explanation is that whenever they initialize their experiment, multiple 
particle pairs are created in some intermediate region. Suppose for simplicity that the par- 
ticles of each pair are identical and have only two levels (Figured]). That is, whenever some 
particle impinges on a detector, it will release a signal of strength A2 (Ai) if the particle 
happened to be on state 2 (1)0; X and Y are just the sum of all such individual signals. 
We will call this theory the microscopic dichotomic model. Alice and Bob could propose 
this model, for instance, to explain why the outcomes of homodyne measurements in usual 
quantum optical experiments follow gaussian distributions. 

We will next show that, in certain circumstances, the microscopic dichotomic model can 
be experimentally refuted, even when no assumptions are made on the values of Ai, A2, the 
number of particles or the number of different species of particles. 

In microscopic dichotomic models, the values of the covariance matrix 7 are given by 

« E^^- E E AiAi,{y;,(a)5(a,a') -p^^(aH(a')}, 

j a=l,2 a'=l,2 

^yy KWBmb,b')-p^{b)p^j,{b')}, 
j 6=1,2 6'=1,2 

7xy oc 5^ AT,. KWAsia, b) - y^(aH(6)}, (2) 

j a,b=l,2 

where p-'^^{a, h) denotes the probability that in one of the pairs of type j, Alice's particle is in 
state a and Bob's, in state 6; p-^(a),p-^(6) represent the corresponding marginal probabilities 
and Nj is the number of pairs of particles of species j. 
On the other hand, the displacement vector satisfies 

(X) cx ^ N,p>^{a)K,. (Y) (X J2 NM<^)Kr (3) 
J i 

^This is the case, for instance, when the particles happen to be photons and the detector is composed of 
a verticaUy disposed polarizer followed by a photocounter. In that case, photons with vertical polarization 
(state 2) will reach the photocounter and thus produce A2 = 1 electrons, whereas no electrons (Ai = 0) will 
be released if the photon was horizontally polarized (state 1). 
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Figure 1: Microscopic model of a macroscopic experiment. N particle pairs travel 
to Alice's and Bob's lab and interact linearly with their detectors. The macroscopic vari- 
ables X, Y observed by Alice and Bob are the coarse-grained sum of all such microscopic 
interactions. 

The values of {X) , (Y) do not carry any information about the original microscopic 
distribution, since we can set them to any value while leaving 7 invariant. Indeed, suppose 
that we add to the ensemble a species of particles such that Ai = 0, A2 7^ 0. Then, setting 
p{l,2) = 1 and varying the value of A2 (or, equivalently, the number of particles of such 
species present in the ensemble), we can set the value of {X) without affecting {X) or 7. 
Likewise, we can set the value of (Y). 

The only information available about the origin of X and Y is thus hidden in 7. At this 
respect, one can check that the following proposition holds: 

Proposition 1. Let X,Y be gaussian variables. Then, (AT, F) can be generated by an en- 
semble of independent pairs of two-level same-type particles iff 



Proof. Let us first prove the necessity of eqs. (jlj). By convexity, it is enough to see that 
the inequalities hold for each of the terms j of the sums in eq. ([2]). Then, written in terms 
of 1), the contribution of each term to 7xx and 7xy (omitting j) is 

equal to (Ai - A2)^Pa(1)[1 -Pa(1)] and (Ai - \2Y[pab{^-, 1) -Pa(1)pb(1)]- Subtracting both 
expressions -and ignoring the non- negative (Ai — A2)^ factor-, we end up with 



Ixx - \1xy\ > 0- 

lYY - \1Xy\ > 0. 



(4) 
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Pa{1)[1-Va{1)] ~Pab{1. 1)+Va{1)vb{1) > 
> pab{1, 1) [1 - Pa{1)] - Pab{1, 1) + Pa{1)pb{1) = 
= Pa{1)[pb{1)-Pab{1,1)]>0. 



(5) 



We have thus demonstrated that 7xx — Ixy > 0. The rest of the inequahties are proven 
analogously. 

To show the opposite implication, suppose that 7 is such that the inequalities (jl]) are 
satisfied, and assume that •jxy > (the case •yxY < is very similar). Then, 

f 1 1 \ , / IXX - IXY \ , . 

7 = 0xk(^1 lj + ( 7VV-7XV j- 

The first covariance matrix on the right-hand side of the above equation can be attained 
by a species of particles following the microscopic distribution pa{c) = Pb{c) = pyis(c, c) = 
1/2, for c = 0, 1. The diagonal matrix corresponds to a situation where the particles compos- 
ing the pair are independently distributed. Adding to the previous ensemble a new species 
of particle pairs with 

(A'l - KTPAim - Pa(0)) = IXX - IXY, 
Ps(0)(l -P'b(0))(A; - X',)' = ^YY - IXY, 
P'ABiO,0)=p'A{0)p'j,{0), 

we thus manage to reproduce 7. 

□ 

Note that, for any r, e > 0, any gaussian bivariate distribution with covariance matrix of 
the form 

Ixx = l/r,jxY = l,lYY = r + e (7) 

violates the inequalities (jlj) as long as r > 1. It follows that not all bivariate gaussian 
distributions can be generated by sums of independent pairs of identical two- valued particles. 
The analysis of macroscopic continuous data can thus give us important clues about the 
microscopic structure that lies underneath. 

On a foundational side, this result implies that Alice and Bob can design a quantum optics 
experiment in order to disprove that entangled gaussian states of light follow a microscopic 
dichotomic model. We will come back to this topic in Section 13. 2[ 

2.1 An application: two-dimensional Brownian motion 

Due to molecular collisions, a particle fioating in a two-dimensional fiuid will experience 
random kicks that will make it move in unpredictable ways. This phenomenon is known as 
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Figure 2: Brownian motion. The dynamics of a particle subject to brownian motion can 
be studied by measuring its displacements X, Y along the x, y axes during some time interval 



Brownian motion, in honor of its discoverer, the botanist Robert Brown. The first theoretical 
explanation of this effect is due to Albert Einstein [TT], who also proposed to use it to estimate 
the size of fluid molecules. 

Experimentally, it is observed that the displacements X, Y of the particle along the axes 
X, y during a given time interval T follow a bivariate gaussian distribution. Hence Brownian 
motion is usually modeled by assuming that, at each time interval t, the particle under 
observation moves a unit of length either up or down and a unit of length either left or right, 
according to the probability distribution {Pt{x^^\ y^^^) : x^^\y^^^ E {il}} (see Figure |2]). If 
Pt is periodic in time, with period Tq ^ T, it is easy to check that the variables 



are distributed according to a gaussian probability density. 

To which kind of macroscopic systems can we apply the former model? Call 7 > G M2X2 
the experimental covariance matrix of {X, Y). From our previous discussion, it follows that, if 
we are to describe the observed phenomena with a random walk, 7 must satisfy conditions dl]). 
In other words, it has to fulfill tT['yWij) > 0, for i,j = 1, 2, where Wij = l2+(— l)*cra; + (— l)"'^"^- 

However, that is not enough. We would also want our microscopic model to be macro- 
scopically isotropic. That is, if we choose to describe the displacement of the particle using a 
different pair of orthonormal axes y', still we should be able to find microscopic dichotomic 
distributions Pt{x^^\y^*^) compatible with our macroscopic observations. This implies that 
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tT{0-fO'^Wij) > 0, for all O G 0(2). 

A bit of algebra shows that this last condition is equivalent to: 

det(7) > ^tr2(7). (9) 

And so a microscopic dichotomic model can only account for Brownian motion in macroscopic 
systems that exhibit the former non-trivial uncertainty relation. 



3 Fixed output structure 

Our aim here is to generalize the problem posed in the previous section. That is, given a set of 
gaussian variables (Xi, X2, we consider the problem of extracting information about 
the underlying microscopic model that gave rise to them. Equivalently, we are interested in 
how different restrictions on such a microscopic model translate to the macroscopic scale. 
Suppose, for instance, that we wonder if (Xi, X„) could derive from several independent n- 
tuples of (i-level particles impinging on different detectors. We conjecture that each detector 
k G {1, ...,n} is simple, i.e., it produces a signal of magnitude when a particle on state a 
impinges on it (we will see later that this assumption can be relaxed). We can thus regard 
each particle impinging on detector A; as a (i- valued microscopic variable Xk G {^a}a=i- 

Suppose now that we have partial knowledge of our measurement devices. More con- 
cretely, imagine that we know the sets of values {A^} up to a proportionality constant. Such 
is the case, for instance, when the response function of each detector is proportional to the 
spin of the incident particle, but we ignore the precise value of the coupling constant (see 
section [3^ . In that case, we could take A^ oc (a — d/2). 

For fixed A = {A^ : k = l,...,n;a = l,...,d}, we want to determine which n-variate 
gaussian distributions can be generated by independent tuples of correlated microscopic 
variables with spectrum proportional to A. In the generic case (whenever each set {A^}a has 
positive and negative values), it is easy to see that the displacement vector Vi = {Xi) can be 
varied at will without modifying the nxn covariance matrix = (XiXj) — (Xj) (X^) of the 
macroscopic system. All information about A must then be contained in 7. Our aim, thus, 
is to characterize those covariance matrices that can be expressed as 

L d 
Z=l a,b=l 

with X^ > being proportional to the number of particle tuples following the distribution I. 
Note that, for any S*, T G M, the families of covariance matrices attainable with the sets of 
outcomes {Af } and {S ■ Xf + T} are the same. 

In order to solve the above problem, as well as some others that will appear along the 
article, we have to introduce the concept of naked covariance matrix. 
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Figure 3: Behavior of a naked detector. A naked detector divides the original beam k 
into sub-beams of particles in the same state a, and then measures the intensity of each 
sub-beam. As indicated in the figure, we can simulate a simple detector by postprocessing 
the outputs of a naked detector. 

3.1 The naked covariance matrix 

Imagine an idealized detector that, when impinged by a c?- level particle, sends a signal of 
magnitude 1 to one counter or another depending on the state of such a particle. If we were 
to send a beam of particles of the same species to this detector, the macroscopic current 
measured on each counter (/c, a) would thus indicate the number of particles of the beam k 
in a particular state a. We will call such a device a naked detector. 

Note that, if the values of {Aj^,...,A^} are known, we can simulate the behavior of a 
simple detector by means of a naked detector, see Figure O Indeed, the single current of 
the former on each round of experiments could be determined by summing up the currents 
registered on each arm of the naked detector, each with its appropriate weight. That is, Xk = 

A^X^. If we were able to characterize the cone of all possible covariance matrices F^j* = 
(X^Xj) — {X^){Xj) generated by n-tuples of (i- valued variables impinging on n different 
naked detectors, we could thus describe the set of covariance matrices generated by simple 
detectors. 

Fortunately, such a characterization is possible. The next technical result will be used 
extensively along the rest of the paper. 

Theorem 2. Let F be an nd x nd matrix. For any pair of lists of outcomes c , c? G {1, d}^ , 
consider the microscopic distributions 

^ n 1 " 

p(ai, 02, On) = - n 5(c, a) + - JJ 5(c , a), (11) 

k=l k=l 
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where S{s, s') denotes the Kronecker delta, and call r(c, cf) their associated naked covariance 
matrices. 

Then T is a naked covariance matrix iff it belongs to the cone generated by {r(c, c?) : 
e {l,...,ci}". 

Proof. The left implication is evident. For the opposite one, it is enough to prove that we 
can recover all matrices generated by a single microscopic distribution pi,...,n(ai, ^2, ^n)- 
In that case, all the entries of the naked covariance matrix are of the form 

r^/ ^ (X^Xi) - {Xp{X',)=p,,{a,b)-p,{a)p,{b), (12) 

where, with a slight abuse of notation, by Pjj{a, a') we mean Saa'Pji^)- Note that we can write 
r"j^ as the difference of two probabilities. Indeed, imagine that we extend our system of n d- 
valued variables (ai, 02, a^) to a system of 2n d-valued variables (ai, 02, a^, a„+i, a2n) 
subject to the distribution p{ai, a„, a„+i, a2n) = Pi,...,n{ai, an) ■ Pi,...,n{an+i, a2n)- 
Then, T'^j' = pij{a,b) — pjj+„(a, 6), and the non-linearity on the probability distribution 
is eliminated. Now p(ai, a„, a„_|_i, a2„) is symmetric with respect to the interchange 
Qi tti+n- Consequently, it can be written as a convex sum of distributions of the form 

g(ai, ...,an,a[, ...,aj = ^6{c,a)6{c ,a) + ^6{c,a)6{c ,d), (13) 
where c,(f G {0,1, ...,d— 1}". The contribution of each of them to T is thus proportional to 

qij{a, b) - qi,j+n{a, b) = ^{S{ci, a)6{cj, b) + 5(c-, a)5{dj, b) - 

-(5(Q,a)5(c;,6)-5(c:,a)5(c„6)}. (14) 

On the other hand, 

Vf^{c, c-) = i{5(c„ a)5(c„ b) + 5(c^, a)5(c;, 6)} - 
-i(5(ci,a) + 5(c-,a))(5(cj,6) + 5(Cj-,6)) = \{qi^j{a,b) - qij+nia,b)}. (15) 

r therefore lies inside the cone of {r(c, cf)}. 

□ 

Now, let 7 be a covariance matrix generated by microscopic variables variables with 
known spectrum {Af}. Since •jij = Yltb=iKiK^ab some naked covariance matrix T, 
Theorem [2] tells us that 7 is inside the cone generated by the covariance matrices 

d 

7.,(c, c-') ^ 4 5^ KK^-^iS, ^) = w{c, d)w{c, df, (16) 

a,b=l 
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where w{c,(^) E M" satisfies w{c,(f)i = A*. — AV 

The problem of determining if 7 admits a microscopic model with known outcome struc- 
ture {A*^} can therefore be cast as a linear program. 

An immediate consequence of this observation is that, as long as ci < 00 and {Af } is fixed, 
there are always bivariate gaussian distributions impossible to attain with such a microscopic 
model. Indeed, suppose that we normalize our two macroscopic variables Xi,X2 in such a 
way that their associated 2x2 covariance matrix 7 satisfies the condition 



tr(7) = 1. (17) 

Then, one can consider which values of 712 are attainable for a fixed value 711. For a 
completely general gaussian distribution, the only condition over 7 is non-negativity. This, 
together with the normalization condition, implies that I712I < \/'yu{l — 711): in this nor- 
malized scenario, the set of attainable points (711,712) is thus a circle, see Fig. HI 

However, suppose that Xi,X2 are generated by independent pairs with outcomes {A^ : 
a = 1, c?}, {A^ : a = 1, ...,d}. According to what we have seen, there exist fi^.s' > such 
that 



7 = ^f^s,ff7{c,c) 



where each 7(0, cf) = ww'^, with 




(19) 



It is obvious that tr(7(c, c')) = implies 7(0, c*) = 0, so these pairs of c,c' do not 
contribute to 7. For the rest, define 7(0, c') = 7(0, c')/tr(7(c, cf)), and note that all 7's are 
such that 7ii,7i2 are in the circumference I712I = a/7ii(1 — 711). Also, we have that 

7 = J]/ic,?7(c,c^), (20) 

and so Ylsc' /^c,? = 1- It follows that 7 is a convex combination of a finite number of points in 
the circumference of the circle. The set of all covariance matrices arising from the microscopic 
model {Xf} is therefore a polygonal inner approximation to the circle corresponding to 
general gaussian distributions, and so the latter contains points separated from the former 
as long as (i < 00. 



3.2 The uniform case: measuring spins 

Let us study a case of physical significance: Alice and Bob perform a quantum optics exper- 
iment in which homodyne measurements are performed on each side of a bipartite gaussian 
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Figure 4: Generic case. A plot of the accessible regions (711,712) of normalized covariance 
matrices compatible with different restrictions on the microscopic structure that gave rise to 
them. The circle with center in (0.5, 0) and radius 0.5 encompassing them all represents the 
set of general gaussian distributions. As you can see, the extreme points of each of these sets 
lie in the circumference. In this example, the inner polytopes represent gaussian distributions 
generated by microscopic variables with values {1/1^, 1/2^..., for d = 2,3,4,5. 

state of light [§, returning a pair of gaussian variables Xi,X2. Having discarded the micro- 
scopic dichotomic model (see Section[2]), they resort to more complicated discrete microscopic 
models in order to explain their observations. That is how they come up with the micro- 
scopic spin model, where Xi,X2 are postulated to be the result of summing up the spins of 
independent {d — l)/2-spin particle pairs, with A^ ^ 1 (note that for d = 2 we recover 
the microscopic dichotomic model). That is, there is a microscopic explanation involving 
only d dimensional quantum systems (Hilbert spaces). We assume that Alice's and Bob's 
detectors behave identically, i.e., whenever a particle with spin component m passes through 
each detector, the corresponding microscopic signal will be equal to g ■ m, where the value 
of the constant g is unknown. We wonder which gaussian states are compatible with such a 
microscopic model. 

Notice that, in this case, we can take {A^} such that Xa — Xb = a — b. Upon imposing 
the normalization constraint ( |T7|) . as long as 7 7^ 0, the covariance matrix will admit a 
decomposition of the form ( 120|) . with 



^For a basic introduction to Quantum Mechanics and Quantum Optics, we refer to [23]. Let us simply 
comment here that these type of experiments are routine nowadays [5] and that the mathematical model 
behind them involves infinite dimensional Hilbert spaces. 



(ci - c[r 



{Ci- C[){C2- C'2) 



(21) 



7ii = 
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Figure 5: Equally-spaced outcomes. A plot of the accessible points (711,712) of nor- 
malized covariance matrices generated by microscopic variables with values {l,...,d}, for 
d = 2,3, 4, 5. The dashed line represents the inequality fl23|l for d = 3. 

Now, we know that the point (0, 0) of the circle is accessible [take Ci = c[]. Our aim now 
is to determine the closest extreme point 7(0, cf) in the circumference. That is, we want to 
find out 

7ti = min{7(c, c'),, : 7(c, ^)ii 7^ 0}. (22) 

c,<? 

It is not difficult to see that 7*^ = Y^z^iziyr- 

Correspondingly, 7*2 = It follows that (see Fig. [5]) the inequality 

1 712 1 < (t^- 1)711 (23) 

holds for all gaussian distributions arising from pairs of {d — l)/2-spin particles, normalized 
or not. 

On the other hand, note that this inequality is violated by all gaussian distributions with 
covariance matrix of the form ([7]) as long as r > (rf — 1). 

From a foundational perspective, eq. fl23|) implies that, for any value of d, there exists 
a quantum optics experiment that proves that gaussian states of light do not follow the 
microscopic {d~ l)/2-spin model. Due to the limitations of current technology, though, this 
claim can only be checked up to a finite value dt- These tests of (in)finite dimensionality in 
quantum optics experiments are rather fashionable nowadays due to their connection with 
security proofs of Quantum Key Distribution systems [7]. 
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3.3 Incomplete output bases 



As d grows, inequality ( l23l) becomes more and more irrelevant. And, actually, one can prove 
that in the limit d ^ oo equally spaced ci-level systems (with positive and negative values) 
can reproduce any multivariate gaussian distribution. Indeed, suppose that we want to ap- 
proximate the gaussian distribution p{Xi, Xn)dX, with covariance matrix 7 and displace- 
ment vector V, by means of d-valued microscopic variables. From previous considerations, 
we only have to worry about reproducing the covariance matrix 7 via microscopic ensembles. 
Now, 7 is positive semidefinite, and so •jki = ■ w\ for some vectors {w^}^^^ C R" such 
that llw'^lloo < L for some L. Define then the vectors {v!'} in such a way that = ^ Lx'^i^J ■ 
It is clear that the new covariance matrix 7^/ = ■ vJ- can be generated by ensembles of 
microscopic systems with outcomes of the form {A^ (x a — d : a = 0, ...,2(i}. Moreover, 
lim^^oo 7 = 7- 

For arbitrary sequences of outcomes {A^ : a G N}, the same argument can be used to 
prove that, as long as 

lim ^^±1^ ^ 1, (24) 

any gaussian multivariate distribution can be approximated by sums of microscopic ensem- 
bles. This includes, in particular, the case Aq = /(a), with / being a non-trivial (i.e., 
non-constant) rational function of a. 

In view of this, it is natural to ask whether there exist sequences of infinitely many differ- 
ent outcomes {A^ : a G N} which do not allow to approximate certain gaussian distributions. 

At first glance, one would be tempted to answer this question in the negative: an infinity 
of possible outcomes gives us infinitely many degrees of freedom to play with. In such 
circumstances, it is difficult to see how, for a given covariance matrix 7, there could not be a 
combination of infinitely many weights that reproduces or approximates 7 at the microscopic 
scale. 

And yet, such is the shown in Fig. [6l when the outcomes happen to be the terms 

of a geometric sequence. The next proposition states it more clearly: 

Proposition 3. Let be the set of all normalized bivariate covariance matrices 7 arising 

from sums of pairs of variables with spectrum {Xa = /i*^ : a G N}. If > ^^/^ , then the 
angle 9 of the extreme points with respect to the horizontal axis satisfies 

Note that, in the limit /i — )• 00, sin(^) ^ (0, 1), i.e., the convex body collapses to the 'bivalued' 
square given by eqs. dl]). Notice also that this results holds too when the outcomes are taken 
from the bigger set {/i° : a G Z}. 

Proof. Any bivariate covariance matrix 7 G is a convex combination of normalized ma- 
trices 7' of the form 
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Figure 6: Geometric outcomes. A plot of the accessible points (711,712) of normalized 
covariance matrices generated by microscopic variables with values {fi°',a G Z}, for /i = 
5, 9, 00. The dashed line delimits the circle of all normalized gaussian distributions. The 
first picture also depicts an extreme point of the circle with sin(6') G 6^. 



7^1 cx ifi^ - fi^y, 7^2 oc ifi' - 7I2 oc (/x" - - (26) 

with a, a', b,b' E N and + 722 7^ 0. We will next prove that, for any such matrix ^ 

(1/(;U — 1)^, (/i — l)^/yU^): the statement of the proposition will then follow from trivial 
calculations. 

We will use reductio ad ahsurdum. Suppose, indeed, that there exist a, a', 6, 6' G N such 
that ^ G (l/(/i — 1)^, (/i — l)^//i^). W.l.o.g. we can assume that a > a',b > b', and then 

the constraint translates as (fi"' — fi'^')/ {fi'' — fi''') G (l/(yU — 1), (/i — l)//u). This automatically 
implies that a ^ a' ,b ^ b' . There are two possibilities: either a = b + A, with A > 0, or 
a = b + A', with A' > 0. 

Let us examine the first one: if a = 6 + A, for A G M, then 

We thus have that a ^ b + A. 

Let us then explore the other option, namely, that a = b + A', with A' > 0. Similarly, 
we have that 

,a ,,a' 1-1 



where the last inequality is due to the fact that A' > 1. 

In both cases we get a contradiction, and so ^ ^ (l/(/^ ~ l)^? (yU — l)^//i^). 

T22 



□ 
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We will say that the set of values {Aq = : a e N} constitutes an incomplete output basis, 
since ensembles of microscopic systems with such an spectrum cannot be used to generate 
arbitrary covariance matrices. 

So far we have been considering symmetric scenarios where = , for k ^ k', i.e., 
where all microscopic systems have the same spectrum. However, the asymmetric case is also 
worth studying, and could play a role in future experiments to reveal hidden quantization. 
For one, in asymmetric scenarios it is not necessary to resort to geometric sequences in order 
to find instances of incomplete output bases of infinite cardinality. Take, for example, the 
pair of spectra {A^ = a : a = 1,2, ...}, {A^ = 1/6 : 6 = 1,2, ...}. Proving that the condition 
hi 2 1 < 7ii holds in this case is left to the reader as an exercise. 

3.4 Imperfect detectors 

All previous results hold in case our detector acts in a deterministic way, that is, when 
our detector assigns a given signal Xa to each particle in state a. However, such is a quite 
unrealistic model for a typical detector behavior. One should rather expect that, depending 
on the state a of the incident particle, our detector sends a signal of strength A, following a 
certain probability distribution Pa(A). Coming back to the case where our detector consists 
of a polarizer connected to a photocounter, the corresponding probability distributions would 
be pi(A) = 5(A), p2(A) = r]5{X — l) + {l—7])S{X), where rj is the efficiency of the photocounter. 

In this section we will show, though, that, at the level of covariance matrices, any proba- 
bilistic detector can be simulated by a deterministic one. Therefore, restricting our analyses 
to this latter case is more than justified. 

Suppose, indeed, that each particle or microscopic variable Xk is evaluated by an imperfect 
detector k, via the following correspondence a — > Pa(A). Let p{Xi, ...,Xk) be the gaussian 
distribution with covariance matrix 7 generated by the microscopic distribution p{xi, 
Defining (/(A))^ = J dXp^{X)f{X) for any function / : a — >■ R, it is straightforward that 

7ij « J2 \PiM - Pi{(^)Pj{b)] (A)L(A)^, if i ^ j, 

a,b 

J2 - PM)PiW)] WiWi' + Ki, otherwise, (29) 

a,a' 

where 

K,^Y.P^{a){{Xy^-{{X)lf)>Q. (30) 

a 

It is therefore clear that if we replace each detector k by the deterministic detector 
a — > (A)^, we obtain a covariance matrix 7' such that 7 — 7' is a diagonal positive semidefinite 
matrix. As long as for each k there exists a{k),h{k) e {1, ...,d} such that (A)^^^) 7^ {^)b{k) 
(and this condition is necessary for the macroscopic variable to be correlated to any 
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other), one can then recover 7 by adding to the previous ensemble a certain number of 
n-tuples of microscopic variables independently distributed with Xk G {{^)^a{k)^ (•^)6(fc)}- 

An immediate consequence of this correspondence between deterministic and probabilistic 
detectors is that relations (jl]) also hold for the latter class, provided that the two detectors 
involved in the experiment are identical. 

3.5 Variable number of variables 

Till now, we have been assuming that our macroscopic variables (Xi, ...,X„) are the result 
of adding up a given number of microscopic variables. Such an assumption, though, may 
not hold in reality: consider, for instance, an experimental situation where the process 
that triggers the production of particle n-tuples is the decay of a microscopic system (say, 
an atom). Imagine also that each atom in our sample decays during the experiment with 
probability v < 1. In that case, our macroscopic signal would show no contributions from 
a particular atom with probability 1 — u. Although we could model the above situation by 
extending the detector response function (i.e., {Xa}i=i — )■ {Xa}i=i U {0}), there is a better 
way to study this scenario. 

Suppose that a given tuple of microscopic variables Xi,...,Xn € {A^} x ... x {A"}, with 
A^ 7^ for any k, a is produced with probability z/. Then one can check that their contribution 
to the macroscopic covariance matrix is of the form 



When z/ = 1, we recover the usual formula, while for z/ ^ 1 we obtain a term proportional 
to (xiXj) = '^abPiji^y^)'^a-Xb- Clearly, an intermediate situation z/ ~ 1 is thus given by a 
conical combination of these two extreme cases. 

If the value of z/ is completely unknown and is allowed to vary between different tuples, 
any accessible covariance matrix 7 will hence admit a conical decomposition of the form 



The problem of deciding if 7 can be generated by ensembles of microscopic variables which 
appear probabilistically in the sum defining the macroscopic variables can thus again be 
formulated as a linear program. Notice that the first (second) summand in fl52]) shall be 
neglected if we postulate that z/ ^ 1 (z/ ^ 1). Unless otherwise specified, along the rest of 
the paper we will always work under the assumption that u = 1. 

4 Free output structure 

We have just studied the case where our macroscopic variables Xi,...,X„ are generated 
by sums of independent n-tuples of microscopic variables with a known structure of values 
{A„ : a = l,...,d,k = l,...,n}. Similarly, one could envision a related scenario where 




(31) 



lij = /^s',c'r"/(c, c^)A^A^^ + nsS{a, Ci)6{b, Cj)\lXl. 



(32) 
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such microscopic variables have d levels, but the concrete values associated to each level 
are unknown, and may even vary between the different tuples. The problem we propose is 
thus to characterize which gaussian distributions can arise from generic averages of d-level 
microscopic systems. 

As it turns out, if we allow the set of outcomes {A^'' : a = 1, ...,d,k = 1, n} for each n- 
tuple I to be completely arbitrary, then any ra-variate gaussian distribution can be generated 
with 2-valued microscopic systems. Indeed, let p{Xi, be a gaussian distribution with 

displacement vector v = (X) and covariance matrix 7. Since 7 > 0, it admits a Gram 
decomposition, i.e., there exists a set of vectors {u^}^^^ C such that 'jij = ■ uK Now, 
define 7^^^- = u\u\, and note that 

7^ 0^ r^''(l,2)Aj,''Af, (33) 

a,6=l,2 

with the values AJ^'^ = , Ag'' = 0. 

7 is thus a conical combination of covariance matrices {7'}[L^ of 2-valued microscopic 
variables. On the other hand, the displacement vector v can be modified by an arbitrary 
amount without altering 7 by adding a constant n-tuple (Afi, ...Af„) to the ensemble. 
The initial gaussian distribution can thus be completely recovered. 

Nevertheless, between a complete knowledge of the outcome structure and a complete 
ignorance there exist natural intermediate situations where the problem of characterizing 
the resulting gaussian distributions becomes non-trivial. Suppose, for instance, that we 
impose the additional hypothesis that the microscopic variables in each n-tuple have the 
same structure, i.e., we postulate that A^'' = A^ for all fc, k' G {1, ...,n} and all /. This is a 
natural assumption when the sequence Xi,X2, represents measurements of the same 

macroscopic variable X at different times, i.e., X^ = Ylh=i ■ In such a situation, it 

is reasonable to postulate that the probability distribution of each microscopic variable x*^'-* 
evolves with time, while its set of possible values remains constant. This scenario may appeal 
to those interested in computational biology, where simple (i.e., not memory consuming) 
discrete mathematical idealizations of biological entities are sought to fit macroscopic time 
series. Also, a closely related assumption will be used in Section [6] to rule out the existence 
of finitely-valued classical models for certain families of quantum experiments. 

From a mathematical point of view, the hypothesis of identical outcomes makes the prob- 
lem non-trivial again: indeed, using Gram decomposition arguments, it is easy to infer that 
any multivariate gaussian distribution p{Xi, ...,X„) can be generated by summing indepen- 
dent n-tuples of n + l-valued identical microscopic systems. The next proposition shows that 
the upper bound + 1 on the minimal number of outcomes is actually optimal. 

Proposition 4. Let p{Xi, Xn) be an arbitrary gaussian distribution with covariance 
matrix 7 generated by convolutions of microscopic distributions p{xi \ ...,Xn^), with x^*^ G 
{X^a^Va^i^k. Then, 
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Note that the former inequahty is violated by general gaussian distributions (take, for in- 
stance, = 2*+-'). 

To prove the proposition, we will need the following result. 

Lemma 5. Let D = {d^}'^^-^ cW^ be a set of vectors of the form = 5{i, Ck) — dj^), for 
c, cf G {1, 2, n}". If the vectors in D are linearly dependent, then there exists a non-null 
vector X G { — 1, 0, 1}*" such that Xkd^ = 0. 

Proof. First, some notation: given a non-null vector d G of the form di = 6{i, c) — 6{i, c'), 
we will call c and c' the occupied indices of d. Also, we will denote by {\i) : z = 1, ...,n} the 
canonical basis of M'^, i.e., |1) = (1, 0, 0, ...), |2) = (0, 1, 0, ...), etc. 

By hypothesis, there exist coefficients {/ifc} such that "^j^fikd^ = 0. Consider the set of 
vectors E = {d^ : /i^ ^ 0}. If there exists e G E such that e = 0, then we have that 1 ■ e = 0, 
and we have finished. Suppose, on the contrary, that none of the elements of E is null. Define 
Xi = 1 and choose an arbitrary element of E (denote it e^), with = — \i2). Since ^ 0, 
there must exist another non-null vector sharing an occupied index with (for otherwise 
'^kf^kd^ ~ would not hold). Call {ii,i2) and (12,13) the occupied indices of e^. Take 
X2 = ef^. If ii = ^3, then J2i=i2^i^ ~ Ni) ~ N2) + K2) — \h) = 0, and we have finished. If 
such is not the case, then there is a third vector whose occupied indices are [i^, 14,). Again 
set X3 = ef^. If ^4 = ii or i^ = 12, then either Y^i=i = or Yl^=2-'^i'^ ~ respectively. 
If ^4 does not equal any of the previous indices, then there exists a different vector e*^ G £" 
with occupied indices (^14^,1^), etc. Since the set of rows is finite, if we iterate the procedure 
at some point we will find a vector whose occupied indices [ik, ik+i) are such that ik+i = ik', 
for some k' < k, in which case we know that Yl^=k' ~ Yl^j=k' \h) ~ Kj+i) =0- ^ 

Proof of Proposition^ In order to prove relation flMl) . it is enough to show that it holds for 
covariance matrices of the form = wuF with the vector w G given by 

Wi = (Ae, - \c{), (35) 

for some c,cf& {1, ra}". Our next step will be to prove that, for any choice of c, cf, there 
exist a pair of disjoint sets A,B & {1,2, ...,n} (one of which may be actually empty) such 
that the non-null vector Si = XA{i) — Xsii) satisfies s-w = Q. 

Indeed, note that w = CX, where Cik = 6{ci,k) — 6{Ci,k). Now, the rows of C must 
be linearly dependent, since there are n of them and all are perpendicular to the vector 
(1, 1, 1). Applying Lemma to C's rows we thus have that there exist a non-null vector 
s G {-1, 0, 1}" such that = 0. 

Now, B{w) can be expressed as 

B{w)= (l--^]\\wf--^{w-vf>0, (36) 
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where v corresponds to the vector Vi = . Let 5*6 { — 1, 0, 1}" be the non-null vector such 
that s'-w = 0, and call A and B the sets of indices where its entries are 1 or —1, respectively. 
Then we have that • s| = | ^jg^2* — ^^£^2*1 > 1, since ^jg^2* and ^jg^ 2* correspond 
to the binary expansion of two different natural numbers. Call s the normalization of the 
vector s. Then, v = fis + y^jj^jp^^T^s-'- , where S"*- is a unitary vector orthogonal to s and 
l/il = \v-s\ > , ^ > ^. 

Finally, we have that 



\w ■ v\ = a/||iTP — /i^ls"*" ■ w\ < y \\v\ 
Substituting in (15^ . we arrive at (151|) . 

□ 



n 



(37) 



4.1 How different are general and finitely generated covariance 
matrices? 

Call Q!^ the cone of all covariance matrices generated by n-tuples of identical d-valued micro- 
scopic variables. Given that both and the set of n-variate general covariance matrices 
are cones, estimating their difference does not make sense unless a scale is fixed. A natural 
way to do so is via normalization, i.e., by dilating each covariance matrix 7 until it satisfies 
tr(7) = 1. Intuitively, this normalization constraint fixes the total amount of noise in the 
system to be equal to unity. 

The next step is to define a suitable distance to quantify the difference of the sets Q^, Q!^. 
One possibility is to use a witness G G Mnxn- That way, for fixed G, a measure of the 
difference between Q'^, would be given by 

r^,(^. ^ max{|tr(Gy)|:yeetr(y) = l} 

^ ^ max{|tr(G7)|:7G^„°°,tr(7) = l}' ^ ^ 

Clearly, K'^{G) < 1, and the lower its value, the greater the distance between the two sets. 
One way to interpret eq. f lM|) is that there exists a witness G G Mnxn such that 

K'\G) = 1 - — ^ -. (39) 

^ ' n(4" - 1) ^ ' 

Indeed, take G,j = 2'+^-^^. G itself is a rank-1 normalized positive semidefinite matrix, 
and so a normalized element of It follows that, for all normalized 7, tr(G'7) < 1, with 
equality for 7 = G. On the other hand, according to eq. (IMI) . tr(G'7') < 1 — „(^Ji_i^ , for all 
normalized 7' e Q^. (^^^ it is not difficult to see that this inequality can be saturated). 

Relation ( 139|) suggests that, although different, Q!^ and Q!^ are exponentially close, and 
thus they would be very difficult to distinguish in practice. However, since general n-variate 
gaussians can be generated by n + 1-dimensional systems, one could argue that the proximity 
between is just a restatement of the fact that n can be approximated by n + 1 in the 
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limit of large n. And actually, if free outcome structure models are to be of any practical 
use, the relevant question is instead how close and are in the limit 1. 

To shed light on this matter, we will follow the lines of [B] and consider an idealized 
scenario where the macroscopic variables form a continuum, i.e., our gaussian variables are 
{Xt : t G [0, 1]}. In these conditions, the covariance matrix of our system shall be replaced 
by a positive semidefinite kernel of the form 7(t, m), such that 

{{XtX^)-{Xt){X^))dtdu= [ ^it,u)dtdu, (40) 

Jn 

where TZ is an arbitrary region of [0, 1]^. We will further assume that our macroscopic 
variables are normalized so that 

1 

-f{t,t)dt = l, (41) 







this condition being the continuum counterpart of demanding that the trace of the covariance 
matrix is equal to 1. 

In this setting, any normalized positive-semidefinite symmetric kernel qualifies as the 
covariance kernel of general macroscopic variables. On the other hand, the kernels generated 
by microscopic variables with spectrum of cardinality d belong to the set 



= cone L{t,u) : ^{t,u) = ^ A,A,(x^(t) - xf (t))(x, 
I k,j=i 



f{u)-xf{n))}, (42) 



where xt (xf ) denotes the characteristic function of the set Ak C [0, 1] {Bk C [0, 1]). The 
sets {Ak}f.^i {{Bk}i=i) are mutually disjoint, and satisfy Uf^iAk = [0, 1] {uf^iBk = [0, 1]). 

In this continuum scenario, witnesses G G Mnxn shall be replaced by bounded symmetric 
kernels G{u, t), and the action of a witness G over an element 7 of will be given by 



G(7) = / G{t,u)-f{u,t)dudt. (43) 

The next result lower bounds the distance between the sets Q"^ and via a simple 
witness. 

Proposition 6. Consider the symmetric kernel G(t,u) = 3tu. Then, 

^ ma.{|GW|:.egJ„S(M^-l} ^ ^ _ 1 

max{|G(T')l :y ea'',/(,y(<,<)*= 1} " - 1) + 2)2 

This proposition has to be understood as a proof that, in the limit n ^ d ^ 1, the 
difference between and decreases at least as the inverse of a quartic polynomial in d, 
i.e., it is not exponentially small. Free outcome structure microscopic models are therefore 
potentially practical. 
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Proof. In quantum mechanics terminology, G can be seen as a normalized pure state G = 
with lip) = / dtt\t). Likewise, any normalized positive semidefinite kernel 'y(t,u) 
can be regarded as a normalized a quantum state 7, and so 

G(7)=tr(G7). (45) 

If we are optimizing the above quantity over general quantum states 7, it is clear that the 
maximum will be attained by taking 7 = G, in which case G{'j) = 1. 

Let us now consider the optimization over Q'^. From equation ( H2l) . it is clear that 
normalized elements of Q'^ are convex combinations of normalized pure states of the 

form 

10)= f dtj2Mxtit)-x'imt)- (46) 
"^0 k=i 

In turn, such states belong to the bigger set 

^ d(d-l)/2+l 

^^{<P:\4>)= dt Yl A.(x^(t)-xf(t))|t),(0|0) = l}, (47) 

where > 0, and, this time, all the sets {Ak, Bj} are mutually disjoint, with Uk{AkUBk) = 
[0,1]. 

An upper bound on the maximal value of ^(7') over the set Q'^ can hence be obtained 
by solving the optimization problem: 



maxtr(7|0)(0|) = max 3 



dt(j){t)t 



(4J 



Now, since t > in the interval [0, 1], we can take all the sets {B^} to be empty. The problem 
we aim at solving is thus equivalent to finding the best approximation of the function t by 
a linear combination of d = d{d — l)/2 + 1 characteristic functions of mutually disjoint sets. 
Since the function t is strictly increasing, it is easy to see that the optimal sets must be 
of the form A^ = [flfe, ak+i], where ai = and a^j^^ = 1. 

Given the intervals [cifc; c^fc+i]? the problem of optimizing {Xk} can be solved by Fourier 
analysis. Indeed, it is clear that, for any vector w and any linear subspace V, max^^v {wl"^) = 
a/ {w\Pv\w), where Py is the projector onto the subspace V. Since, by definition of the sets 
{Ak}, the functions {Xfc (^)}fc=i cire orthogonal, we have that 



{ilj\P^\^) = J24i(^k+i-ak), (49) 

where 



k=l 
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The problem thus reduces to optimize the intervals [a^, a^+i] so as to maximize 3/(a), where 



Finally, note that 



/(«) = E 



A;=l 



(51) 



fc=lV'^fc+l '^feJ 



/(a) 



E 



fc=i 



"1+1-4 

3{afe+i-afe) 



Ed 
k=l 
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(52) 



Renaming = a^+i — at, our aim is now to minimize \\v\\3 over all vectors w G M'^ such that 
llvlli = 1. This is a well-known problem in Banach space theory: the optimal vector must 



satisfy \vk\ = j. Hence we end up with: 

max{tr(77') : 7' G Q^, f dt^'{t,t) = 1} < min(^|P^ 

Jo A 

Substituting d = d{d — l)/2 + 1, we arrive at 



1 - 



4(i2- 



(53) 



□ 



Note that we can reinterpret eq. ( 14^ as a lower bound on the trace distance between a 
certain normalized element of and the set of normalized elements of Q'^. Indeed, as we 
saw during the proof of proposition [6l G G and G(t,t)t = 1. On the other hand, for 
any '~f & Q'^ satisfying the normalization condition, we have that 



|G'-7||i>2[l-G'(7)]> 



(54) 



Kd- 1) + 2]2' 

where the first inequality derives from the relation 11^ — 7111 > trKG — 7)(2G' — I)}, and the 
second inequality follows from eq. (jB]). 



4.2 General algorithm 

Now that we have proven that the general problem makes sense, the next question to ask 
ourselves is how to solve it in general. That is, given an n-variate gaussian probability 
distribution with covariance matrix 7, how can we determine if it can be generated through 
independent ra-tuples of microscopic variables {x[,X2, ...,x^j^) G {A^^, A^}*^? 

A blind application of Theorem [2] leads one to the following semidefinite program |22j : 



min 

s.t. 

^^^• = EEr(c,C-')^-^a6(c,C-'), 

c,<? a,b 

Z(c,c')>0. (55) 
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Indeed, one can see that the covariance matrix 7 generated by an arbitrary microscopic 
distribution p{ai, 02, ...) is equal to 

= Ea,b Ec,^ f^c,^^t^ic, ^)KXb = 
= Ec>- Ea,b^f!ic,^)Zal>{S,^), (56) 

with Zab{c,c') = fJ'c,c'^a^b- By convexity, we thus have that program fl55l) defines a set of 
matrices that contains Q!^. 

Conversely, let 7 = ^ F^j* ■ Zab, where F is a naked covariance matrix and Z > 0. Then 

one can Gram-decompose Z as Zab = Yli=i ^aK^ so 7 is generated by the d n-tuples 
{x\,...,xl)e{X\,...,Xiir,t = l,...,d. 

It follows that program (155|) completely characterizes the set Q!^. 

Program (!55|) is, though, unnecessarily inefficient. Notice, for instance, that the covari- 
ance matrices do not change if we fix Ai = 0. Consequently, the matrices Z{c, c) can be taken 
of size (d—l) X [d — 1). Notice as well that certain pairs of points (c, &) and {d, d') actually 
codify the same information: for example, the matrices w{c, (f)w{c, cf)^ and w{(f, c)w{(f, c)^ 
are identical. Also, in the case d = n = ?> the vectors 

wi = A2 - Ai, W2 = Ai - A3, W3 = A2 - A3, (57) 

and 

U7i = A2 - Ai, W2 = A3 - A2, W3 = A3 - Ai, (58) 

although seemingly very different, actually contain the same information (namely, that wi 
and W2 are free and wz = wi + W2). Likewise, the matrices generated by a vector of the form 

wi = A2 - Ai, W2 = 0, W3 = A2 - Ai (59) 

are among the ones generated by any of the other two vectors. 

How to avoid this redundancy? The solution is to divide the set of possible points 
{!,..., d}"" X {1, d}"' into different classes and choose just one representative (if any) of each 
class for the semidefinite program. The classification we propose is based on the following 
fact: given c, cf, one can always find a unique parametric representation for w of the form 
w = Pt, with P being an n x [d — 1) matrix with entries in { — 1,0, 1} and satisfying the 
properties: 

1. The first non-zero row of P is (1, 0, 0). 

2. Let N{j) denote the greatest from left to the right non-zero index of rows 1,2, 
Then, N{j + 1) < N{j) + 1. 

3. If N{j + 1) = N{j) + 1, then P.+i,^ = 6k,NU+i). 
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The intuition behind this canonical form is that those rows j where N{j) = N{j — 1) + 1 
are independent parameters, while any other row is a linear combination of the first N{k) 
'independent rows'. From the proof of Lemma it is clear that the entries of P have to be 
O's and ±l's. 

Notice that if for some matrix P, N{n) < d — 1, then the covariance matrices generated 
by P are a subset of those generated by a modified matrix P'{P) where we have made 
d — 1 — N{n) dependent rows independent. In sum, we have only to consider those matrices 
P satisfying conditions 1-3 and generated by a couple of points {c, c'} such that N{n) = d—1. 
We will call the set of all such matrices. Note that V^~^^ = {I}. 

To clarify these ideas, consider the case n = d = 3. Then, one can check that the elements 

'I 



of Vn are: 



Pi=\^l j,Pi=\^^ j,P;^'=\^0 1 j , (60) 
with i,j G {0, ±1}. And, consequently, a 3 x 3 covariance matrix 7 belongs to ^| iff 

i=0,l,-l ij=0,l,-l 




^,3 
3 • 



for some 2x2 positive semidefinite matrices Zi, Z^, Z- 

How does this set of covariance matrices look like? Figure ([7]) shows a plot of the regions 
= {({^7^i,^f7t;2) : tr(7) = 1,7 e Gl}, for d = 2,3,4, with = (1,2,4), = (4,2,1). 
We used the MATLAB package YALMIP [H] in combination with SeDuMi [21j to perform 
the numerical calculations. The sets and Q^, although very different from Q^, seem quite 
similar to each other in this two-dimensional slice. 

Let us finish with a remark on the applicability of the former results. Notice that program 
fl55|) is only valid if the events that generate the tuple production occur with probablity u = 1 
during the course of the experiment. In case u is completely unknown, the expression for 
the covariance matrix 7 shall be complemented with another summand of the form 

Y,S{a,Ci)S{b,Cj)Yabic), (62) 

a,b,c 

with Y{c) > for all c G {1, ...,d}^ (see Section [X^ . Moreover, in the regime z/ ^ 1, only 
this summand should appear in 7's decomposition. As before, by exploiting redundancies, 
one can simphfy eq. (162!) to a great extent. 



4.3 The dual problem 

Given that Q!^ is a cone, an alternative way to characterize it is through linear witnesses, 
i.e., n X n matrices G with the property that tr(G'7) > for all 7 G Q'^, but such that there 
exist general covariance matrices 7' G Q'n'^^ with tr(G'7') < 0. According to what we have 
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Figure 7: Free outcomes. A slice of the accessible regions of normalized 3x3 covariance 
matrices generated by tuples of two, three and four-valued microscopic variables. 

seen, in order to certify that G is indeed a witness for Qf^ all we have to do is to check that 
G ^ and that P'^GP > for all P e The proof is trivial. 

A comparison between the structure of and those of the sets of classical and quan- 
tum correlations is in order. Given a linear functional {Bjl}, in order to verify that 
'Yliabxy-^ai)^i^^^\^-'y) — ^ holds for classical distributions P{a,b\x,y)j one has to calculate 
a finite number of real numbers (the evaluation of B in the extreme points of the classical 
polytope) and check that all of them are greater or equal than 0. On the other hand, to 
prove that the inequality holds true for quantum distributions, one would have to show that 
the (infinitely many) Bell operators that result when we substitute P{a, b\x, y) by projectors 
of the form §>> are positive semidefinite. 

In contrast, certifying that G{j) > for all 76^^ amounts to evaluate a finite number 
of matrices and check that each of them is positive semidefinite. At a purely mathematical 
level, the set is thus very interesting, since its behavior is intermediate between the sets 
of classical and quantum correlations. 

5 Problem complexity 

Despite our simplifications, the computational power required to run the algorithms pre- 
sented above to solve the ^^-membership problem grows very fast with n. This makes one 
wonder whether it is possible to find simpler schemes to characterize finitely-generated gaus- 
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sian distributions. In this section we will show a negative result in that direction: deciding 
whether a given matrix 7 belongs to is a strong NP-hard problem. That is, there exists 
a polynomial q{n) such that the membership problem for the set C = {7 G ^^|tr(7) = 1} 
cannot be decided with precision better than in euclidean norm unless P = NP. 

This result is a trivial consequence of two recent contributions. One, due to Liu [131 
Proposition 2.8], reduces the problem to linear optimization over C with inverse polynomial 
precision. Verifying the hypothesis of Liu's result is trivial in our particular case. Then, one 
can invoke the following result of Bhaskara et al. [21 Appendix B] and notice that QPRATIO 
is nothing but a linear optimization over C. 

Proposition 7. There exist two absolute constants a > P such that, given a matrix A with 
coefficients bounded by 1 and zeros in the diagonal, it is NP-hard to distinguish between a 
solution > a and a solution < 13 in the following problem, called QPRATIO: 

max ^ — K 

{1-1,0}" 

As a final remark, our manuscript shows a natural context where the maximization 
problem QPRATIO appears, which can foster the recently initiated research on this problem 
within the computer science community. 



6 Connection with quantum non-locality 

There exist situations where several macroscopic variables are at stake, but nevertheless the 
probability density p{Xi, . . . , X„, Fi, . . . , F„) is not experimentally accessible (actually, it 
may not even exist). This is the case when two parties, call them Alice and Bob, are allowed 
to interact with their respective ensembles before a measurement is carried out. In such 
scenarios, for each possible pair of interactions z, t, Alice and Bob will be able to estimate 
the marginals p{Xz, Yt\z, t)dXzdYt, see Figure [HI Moreover, if Alice and Bob's operations are 
space-like separated, then ^(X^l^;, t) = p{Xz\z,t') {p{Yt\z,t) = p(Yt\z' ,t)), for all t,t' {z,z'). 

Suppose now that all such probability densities happen to be gaussian. A plausible 
explanation Alice and Bob may come up with is that, whenever they perform an experiment, 
there is an event in some intermediate region between their labs that produces ^ 1 
independent pairs of d-\eveled particles. By conservation of linear momentum, two particle 
beams are thus produced, one directed to Alice and another one to Bob. The action of 
Alice's (Bob's) interaction z (t) affects individually the state of each particle in her (his) 
beam. After such an interaction, the particles impinge on a simple detector, that produces a 
state-dependent signal for each incident particle, and Alice and Bob's readings are precisely 
proportional to such a sum of signals. 

While the particles in each pair will not be necessarily identical, it is not unreasonable to 
assume that Alice's (Bob's) interaction does not change the nature of such particles, but only 
their level probability distribution. That is, even though 7^ A* , it is natural to suppose 
that \l = XI , \l = \l . We will introduce the notation GnAn% denote the set of gaussian 
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Figure 8: Bipartite macroscopic quantum experiment. Alice and Bob receive each a 
beam of correlated particles, with which they are allowed to interact in some ways z and 
t before they reach the detectors. Consequently, they are able to estimate the probability 
densities P{Xz,Yt\z,t)dXzdYt. If their measurements are sufficiently coarse-grained (but 
precise enough to detect fluctuations in the values of Xz,Yt), this set of gaussian marginal 
distributions will always admit a local hidden variable model, that is, a classical explanation. 



distributions with ua {it-b) macroscopic variables associated to Alice (Bob), generated by 
c?A('^s)-identically-valued classical systems. 

Due to Macroscopic Locality (ML) [16j, for any set of pairs of gaussian distributions 

{p{X,,Yt\z,t)dXM},,t (63) 

generated in a quantum bipartite experiment, there is always a global probability distribution 
p{Xi, Xn, Yi, Yn) that admits {p{Xz, Yt\z, t)dXzdYt}z,t as marginals. 

One way to interpret this result is that the outcome of any macroscopic quantum exper- 
iment can be explained by a classical particle model, where each classical particle reaching 
Alice's (Bob's) lab produces a current proportional to X^ (l^t), according to the distribution 
p(Xi, X„, Yi, Yn) El. This classical theory, though, seems unnecessarily convoluted: note 
that, independently of the number of levels dq in the original quantum experiment, the cor- 
responding classical particles have a continuous spectrum. One thus wonders if there exists 
a simpler mapping between quantum and classical macroscopic systems that preserves the 
finite cardinality of the microscopic spectrum. In other words, given a macroscopic quantum 
experiment with dq < oo, is it always possible to find a classical description with dc < ool 

In this section we will answer this question in the negative. More concretely, we will prove 
that, in order to reproduce quantum macroscopic experiments involving many copies of the 
maximally entangled two-qubit state [dq = 2), one needs to conjure up classical models with 
particles of infinitely many levels [dc = oo). So, in a sense, quantum non-locality leaves some 

■^This is no longer true at the microscopic level. The existence of bipartite distributions that can be 
generated by quantum means and not by classical means, called quantum non-locality, was already theoret- 
ically observed by Einstein, Podolsky and Rosen in 1935 [T2] and, based on the results of Bell T, proved 
experimentally by A. Aspect et al. [T]. Since the pioneer work of A. Eckert [TU] , quantum non-locality is 
nowadays the key resource to have device-independent security in quantum cryptographic protocols, such 
as random number generation of quantum key distribution |15j . Other uses of quantum non- locality 
include i.e. exponential communication savings in multipartite computation [5]. 
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signature in the macroscopic level. We refer to [19] for basic quantum information notions 
and notations that will be used in this section. 

Take the following functional, based on the chain inequality 

n n— 1 

= + - (64) 

k=l k=l 

where = — (Xf^), Y^. = Y^. — (Xk)i and consider the problem of minimizing it under the 
assumption that 

n n 

Y.^Xl) + Y.^Y^) = ^- (65) 

k=l k=l 

Then we have the following proposition. 

Proposition 8. For any n, the minimal value of ( [5^ compatible with Macroscopic Locality 
is equal to g^ = — cos (^), and can be attained by performing equatorial measurements over 
the two-qubit singlet state :^(|0)a|1)_b — |1)a|0)_b). 

Proof. From ML we know that, for any set of bipartite marginal distributions, there exists 
a classical model reproducing the observable second momenta of our macroscopic variables. 
It follows that optimizing Gn{{{XjYk)}) against all two-point correlators compatible with 
ML is equivalent to an optimization over all trace-one positive semidefinite matrices 7. By 
convexity, we can thus take 7 = w{w)'^ , with w G M."^^, \\w\\ = 1. 
Consequently, minimizing Gn is equivalent to optimize 

n n—l 

Gn{w^ w^) = wtw^ + ^ - wtw^ (66) 

k=l k=l 

over all vectors w^,w^ G M" such that Hw'^P + H'W^^IP = 1- 
So let us first optimize w^. We have that 



> - + + K - O^ll^^ll- (67) 



Our new task is thus to maximize 



n 

G'{w^) ^ J](^f + wl,r + [wf - w^r, (68) 

k=2 

Note that we can express G'{w^) as {w^)'^G'w^ , where G' is an n x n matrix whose 
non-zero components are 
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6-^ = 2, for^ = l,...,n 

C'm+i = Gi+i^i = 1, for i = 1, ...,n - 1 

G[,n = Gn,i = -1. (69) 

Gathering intuition from circulant matrices, we try the ansatz = e'"^ for the eigen- 
vectors of G'. One finds that 

{G'w^)k = 2e'"'= + e'"^+i + e'^'^"^ = 2[1 + cos(a)]w;f , (70) 
for k = 2, ...,n — 1. However, 



{G'w^)n = 2e'"" + e'"("-^) - e'". (71) 

We circumvent this problem by imposing that e*"'^""'"^) = — e*". This leads to a = (2/ + 
l)7r/n, with / = 0, 1, — 1. This ansatz thus allows us to recover all eigenvectors of G'. 
The greatest eigenvalue is 2(1 + cos(7r/n)), and it is attained by the vectors {wk oc e*'^'^/"} 
and {wk (X e~*'^'^/"}, or, equivalently, by any vector of the form {wk = \/2 sm{kTr / n + 5)}. 
The space generated by such vectors will be denoted by W, and will play an important role 
in the proof of Proposition [TOl 

Coming back to our present problem, we have that the minimum value of flM|) upon no 
restrictions on the gaussian distribution is equal to 



1 + cos ( — ) max ||tZ;'^|| = — cos ( — ). (72) 

\nJ i \\wA\\2+\\^B\\2^i \2nJ 

To complete the proof we thus have to show that the above value can be attained in a 
macroscopic experiment where many particle pairs in the singlet state are produced and Alice 
and Bob's interactions correspond to projective equatorial measurements of such particles. 

It is straightforward that the measurements 



Xk = cos(V'fc)cr^ + sm{iJk)(^z, Vk = cos{4>k)(T^ + sin(0fc)o"^, (73) 

with ijjk = ■^^^r^! 0A; = ^^^2n^'''^ ^11 both sides of a singlet state saturate the bound (7*, provided 
that we assign the values {^5 — ^Sich possible outcome. 

□ 

We just saw that bi-valued quantum variables suffice to attain the minimum value (7* 
of Gn{{{XjYk)}). How does this number relate to the minimal value attainable through 
classical systems of dimensions dA,dB'? The following proposition studies the limits of the 
(Ia = 00, ds = 2-scenario. 
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Proposition 9. The minimal value of {tr{Gn'y) '■ tr{'~f) = 1,7 G G!i^n} ^■^ 




Note that < g^ioi n > 2. Moreover, lim„^oo tt^ = 0. 

Proof. Following the proof of Proposition [HI we just have to maximize (!68|) for G {0, ±A} 
as a function of Note that the minus sign in (168|) between entries 1 and n destroys 

translational invariance. However, we can displace it to any position k of the chain by 
performing the change of coordinates wf — )■ —wf, for i < k. Now, suppose that G 
{±A, 0}" has exactly s > 1 null entries, choose the first site k such that wj^ = and 
7^ 0, and place the minus sign between k and k — 1. That way, for the purposes of 
maximizing G'{w^), we can always choose the non-zero components equal to A. 

Imagine now that all sites except site k are equal to A, giving a value G'{w^) = 4A^(n— 1). 
We are now going to place the remaining s — 1 zeros in the chain starting from site k and in 
increasing order in such a way as to maximize the value of G'{w^). While 'nullifying' each 
site i, we find two different situations: 

1. The sites (i — 1, z, z + 1) are of the form — A — A and we substitute the middle one 
by a 0, i.e., we will have — — A. In that case, the value of G'{w^) will decrease by 
an amount of 4A^. 

2. The sites (i — 1, i, z + 1) are of the form A — A — A, we substitute the middle one by 
a 0. That way, we arrive at A — — A, and the corresponding decrease in G'{w^) is 
equal to 6A^. 

It follows that, in order to maximize G'{A) for a fixed number of zeros s, with n > s > 1, 
the best strategy is to place them all one after the other. The maximal value of G'{w^) is thus 
(4- The best situation is hence s = 1, in which case G'{w^) = (4 - 

It only rests to consider the case s = 0. But it is immediate that then the maximum value 
of G'{w^) equals the slightly smaller (for n > 2) value (4 — 8/n)||w^|p. 

The smallest value of (IMI) under the assumption that Bob's microscopic variables are 
classical and bivalued is therefore 
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min - /maxG"(M)||w^||||w-^|| = -Jl - — ^ — -. (75) 

H|2+||^s||2=i Y,|.s,|=i ^ III! II Y 2(n-l) ^ ' 

□ 

The case = oo, > 2 is much more difficult to analyze, and so this time we can 
only provide gross lower bounds h'^ for = min{tr(G„7) : tr(7) = 1,7 G Q'^n}- Such 
bounds, though, satisfy /ij^ > (?* for d ~ ©(-y/n), thereby proving that no finitely- valued 
classical model can account for the macroscopic correlations we observe when we perform 
polarization measurements over different ends of a double beam of many photons in the 
maximally entangled state. 
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Proposition 10. Let \d{d - 1) + 3 < [f ] . Then, > g*^. 

Proof. From the proof of Proposition [HI we have that, in order to attain (7*, Bob's vector 
has to be parallel to a vector in W , the 2-dimensional space of eigenvectors with minimum 
eigenvalue. Hence, if we manage to prove that cannot approximate any such vector, we 
are done. 

Suppose then that — 1) + 3 < [|], and call u the (normalized) projection of 
onto W . Then, the number of different non-negative (or non-positive) entries {Aj — \j : 
i,j = l,...,d} of is at most [|] — 2. Now, any unitary vector m G is of the form 

Uk = y^sin (^fc -|- 5'). Identifying the (n + I)*'' entry with the first entry, we thus have that 

u must have a sequence of consecutive entries of the form ±^ sin (^^j + 6) , for j = 0... [|] —2, 
with < 6 < ^. It follows that there exist at least two entries k, I such that = wf and 

Uk = ±y^sin(7rmfe/n -|- (5), ui = ±^J^sm{Tcmk/n + 5), with 171^,1111 G N, 7^ mi. The 
minimum of the quantity \\w^ — is thus lower bounded by 



> |(m; - UkY > mmmi^jLmi ^ [siu (mfc^ + 6) ~ sin [mi^ + 6)f > 

>^sin2 (It) «in^(^)>o, 



1 , 



(76) 



where the last inequality follows from taking = [|] — 2, m; = [|] — 3 and applying some 
trigonometric identities. 

It follows that ^ W, and so the minimum g^ cannot be attained. It is straightforward 
to carry on this argument in order to find a horribly complicated lower bound on the difference 



□ 



In view of this last witness Gn, whose maximal value was attained by quantum mechanical 
bivalued variables, one may be tempted to think that quantum dichotomic systems suffice to 
reproduce all gaussian marginal distributions arising in a bipartite scenario. However, such 
is not the case. Consider a scenario where Alice has only one interaction to choose from. 
Then we have the following result. 

Proposition 11. Let Alice and Bob's macroscomc observables arise as sums of microscopic 
variables subject to the no-signalling constraint^ only, with dA = 2 and ds = n. Then, the 
macroscopic variables satisfy: 



k=l 



1 

n 



k=l 



> 0. 



(77) 



"'The no-signalling condition states that the marginals cannot depend on the interaction chose by the 
other party. Otherwise, assuming the parties are space-like separated, one would have arbitrarily fast -in 
particular faster the light- communication. Obviously, quantum mechanics is a no-signalling theory. 



31 



This inequahty can be violated by the classical matrix 
with w"^ = 1, = y^^^2''~\ for k = 1, n. 

Proof. Since Alice has only one interaction available, there exists a global probability distri- 
bution for the variables {x,yi, ...,yn) (take, for instance, p(a, 61, 6„) = p(^)n^i YYi=iP(^^^i))- 
It follows that we can regard the scenario as classical and apply the techniques developed so 
far. The problem thus reduces to minimizing the expression 

{v-w'')w'^, (79) 

where Vi = 2^~^,i = l,...,n, for a fixed value of Hu^^P + (w^)^. Following the proof of 
Proposition HJ we have that ± s, with Si = Xcii) — Xoij), for two disjoint sets of 
outcomes C and D such that C U D 7^ 0. Using the same arguments, one concludes that 

\v ■ < y^ji^^|p"-^||w^||, and so we arrive at (1771) . 

□ 



7 Conclusion 

In this paper we have shown that different microscopic spectra lead to different restrictions 
on the macroscopic probability distributions that result out of the convolution of a big num- 
ber of independent microscopic systems. To put it in another way: it is possible to extract 
non-trivial information about the output structure of a microscopic system just by analyzing 
the corresponding gaussian macroscopic variables. We have shown that characterizing the 
set of all gaussian distributions generated by microscopic systems with a given (finite) spec- 
trum can be formulated as a linear program. We have also proven the existence of infinite 
dimensional outcome structures that do not allow to recover the set of all gaussian bivari- 
ate distributions. In the case where the spectrum of the microscopic variables is not fixed 
a priori, we have shown that n + 1 outcomes are enough to recover all n-variate gaussian 
distributions, and this number is optimal. We have demonstrated that the problem of char- 
acterizing gaussian distributions generated by independent tuples of d-valued microscopic 
variables with free outcome structure can be cast as a semidefinite program. Furthermore, 
we studied the complexity class of characterizing the witnesses of the set of gaussians gener- 
ated by dichotomic variables. Despite the fundamental nature of all these problems, to our 
knowledge, this line of research had never been previously considered in Statistics. 

We used the above results to prove that bipartite gaussian states of light cannot be 
regarded as ensembles of pairs of yet-to-be-discovered d-level particles. More concretely, we 
showed that, for any value of there exists a quantum optics experiment that proves that 
gaussian states do not follow the microscopic {d — l)/2-spin model. Likewise, we saw that 
classical models aiming at describing certain macroscopic quantum experiments involving 
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dichotomic quantum systems necessarily involve infinitely-valued classical variables. From a 
foundational point of view, these two results justify the whole study. We believe, though, 
that the range of applicability of the techniques developed in this work does not end here, 
and that they will eventually find a use in other fields. 

From a strictly mathematical point of view, this work leaves open some problems that 
we believe are important to address: 

1. In Section El we showed a witness for Q^^^^'' that was violated by quantum di- 
chotomic systems. However, the exact value = min{tr(G„7) : tr(7) = 1,7 G G^n} 
could not be computed exactly. Is it possible to improve the bounds for Alter- 
natively, is there a better witness to separate the quantum bivalued set from classical 
many-valued models? Note that Gn is not very robust against detector noise. 

2. If very high precision measurements are available, we may obtain (non-zero) estimates 
of higher order cumulants of Xi, X2, X„. Can we use these estimates to extract 
more information about Xi, ...,a;„? 

3. In this work, we have characterized the set of gaussian distributions arising from clas- 
sical finitely-valued systems. Can our results be extended to the quantum case? That 
is, given a set of marginal gaussian distributions P{Xz,Yt\z,t), is it possible to de- 
termine if they can be generated by microscopic quantum systems? In the fixed out- 
come case, it is easy to think of a hierarchy of SDPs to bound such a set of gaussian 
distributions: following the derivation of Theorem [21 consider four-partite quantum 
distributions p{a,b,a' ,b'\z,t, z' ,t') which are invariant with respect to the interchange 
(a, b) -H- (a', b'). Such distributions can be bounded by sequences of SDPs [171 [THl [9], 
and the cone of naked covariance matrices generated by p{a,b\z,t) can in turn be 
bounded by a linear expression on p{a, b, a', b'\z, t, z', t'). Does this approach converge 
to the desired set of gaussian distributions? If not, one could consider extensions of 
p{a,b\z^t) to 2n parties and invoke some (to be discovered) de Finetti theorem for 
quantum correlations. Does this new approach converge? Even better: is it possible 
to characterize 'quantum' covariance matrices with just one SDP? 
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